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A model is presented for the quantum transport of electrons, across finite atomic wire nanojunc- 
tions between electric leads, at zero bias limit. In order to derive the appropriate transmission and 
reflection spectra, familiar in the Landauer-Biittiker formalism, we develop the algebraic phase field 
matching theory (PFMT). In particular, we apply our model calculations to determine the electronic 
conductance for freely suspended monatomic linear sodium wires (MLNaW) between leads of the 
same element, and for the diatomic copper-cobalt wires (DLCuCoW) between copper leads on a 
Cu(lll) substrate. Calculations for the MLNaW system confirm the correctness and functionality 
of our PFMT approach. We present novel transmission spectra for this system, and show that 
its transport properties exhibit the conductance oscillations for the odd- and even-number wires 
in agreement with previously reported first-principle results. The numerical calculations for the 
DLCuCoW wire nanojunctions are motivated by the stability of these systems at low temperatures. 
Our results for the transmission spectra yield for this system, at its Fermi energy, a monotonic 
exponential decay of the conductance with increasing wire length of the Cu-Co pairs. This is a 
cumulative effect which is discussed in detail in the present work, and may prove useful for appli- 
cations in nanocircuits. Furthermore, our PFMT formalism can be considered as a compact and 
efficient tool for the study of the electronic quantum transport for a wide range of nanomaterial wire 
systems. It provides a trade-off in computational efficiency and predictive capability as compared 
to slower first-principle based methods, and has the potential to treat the conductance properties 
of more complex molecular nanojunctions. 
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I. INTRODUCTION 

Current technological needs motivate the intensive re- 
search of the properties of materials at the nanoscale. 
One of the most important areas in this respect at present 
concerns nanoelcctronics The great interest in this 
domain arises from the potential reduction of the size 
of the circuit components, maintaining their quality and 
functionality, and aiming at greater efficiency and storage 
characteristics for physical devices 0- 

Among the various physical components which con- 
stitute the nanoelectronic devices, the nanojunctions 
therein are of crucial importance and will become in- 
creasingly so in the future [H , 3 • Such low-dimensional 
elements, which can act e.g. as switches or gates 0, 
are considered to be key circuit components connecting 
the larger nanostructurcs. At present they are experi- 
mentally prepared and investigated by means of several 
techniques, such as the mechanically controllable break 
junction methods [H, [1], feedback-stabilized break junc- 
tion technique Q , atomic force microscopy [l| , scanning 
tunneling microscopy_^ j j as well as transmission elec- 
tron microscopy Q, 
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Atomic wires represent a particular class of such nano- 
junctions in nanoelectronic circuits. The conductance 
properties of atomic wires are the most interesting fea- 
tures of these systems, they depend on the materials used 
to fabricate them and on their structural properties. For 
example, the conductance of the monatomic wire nano- 
junctions does not stay constant or decrease when the 
length of the wire increases, rather it oscillates as a func- 
tion of the number of the wire atoms. The periodicity of 
the oscillations may vary depending on the type of the 
atomic wire. It has been shown that, in the monovalent 
wire nanojunctions composed of alkali (Na and Cs) or 
noble metals (Ag, An and Cu), the conductance oscilla- 
tions exhibit the two-atom periodicity @ , [l^ , [ll| , [l^l ■ 
In contrast the Al wire nanojunctions present the oscil- 
lations with a period of four or six atoms, depending on 
the considered interatomic distances between consecutive 
sites [l3| . [13. As well as in Al wires, the conductance os- 
cillations have been observed in other multivalence wire 
nanojuctions, as platinum and iridium [l^, carbon [l6| . 
[13, and siHcon systems fisj. figj. 

Understanding the electronic transport properties of 
such nanojunctions is hence of particular importance. 
The electrons which contribute to the transport via nano- 
junctions present characteristic wavelengths comparable 
to the size of these components, leading to quantum co- 
herent effects. The properties of the nanoelectronic de- 
vice and its functionality may be greatly affected or even 
built on such nano junction quantum effects [13], and can- 
not be described in the framework of the classical regime 
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[21| . In particular, the transport characteristics are usu- 
ally described employing the scattering theory for elec- 
tronic excitations in nanostructures |22|- The relation- 
ship between quantum scattering and the electronic con- 
ductance is initially provided by the work of Landauer 
[2^, together with its generalization to the case of mul- 
tiple scattering as developed by Biittiker [2J]. The de- 
termination of the scattering transmission and reflection 
probabilities for the electrons as quantum particles in 
nanostructures is consequently a key concern. 

The scattering probabilities can be determined thanks 
to a number of theoretical methods which have a common 
objective in the framework of the Landauer-Biittiker for- 
malism, and can be employed by integrating in their ap- 
proach, the density functional theor y (D FT) j25j . (26| or 
the tight-binding model (TB) [111, The most pop- 

ular methods arc the matrix Green's function method 
by Ando [2^ or the non-equilibrium Green's function 
formalism (NEGF), also known as the Keldysh formal- 
ism [3^. However, there has been a debate whether 
the Ando's formalism is complete [3l|, and as has been 
pointed out, the combination of the Keldysh formalism 
and DFT is still very often computationally demanding 
[HI J [13: motivating intensive research for solutions to 
this difficulty. An interesting discussion of different types 
of the optimization techniques may be found in [25j , [26j , 

In this paper we present a compact and efficient theo- 
retical method to calculate the electronic quantum con- 
ductance, in the framework of the Landauer-Biittiker for- 
malism, for the class of nanojunctions made up of finite 
atomic wires of conducting atoms. Our model calcula- 
tions are based on the method of the phase field match- 
ing theory (PFMT), which was originally and intensively 
developed for the scattering and transport of phonon and 
magnon excitations in nanostructures [33| . While at first 
sight, the formulation of the PFMT may be expected 
to be formally similar for different types of excitations, 
i.e. replacing here the dynamical equation for phonons or 
magnons by the Schrodinger equation, there are however 
significant technical dissimilarities which play an impor- 
tant role. For example, we can refer to the fact that in the 
electron scattering problem one has to consider allowed 
symmetries of the atomic orbitals, which is not the case 
for phonon and magnon excitations. 

To fully benefit from the PFMT method, our model 
dynamics arc described in the framework of the linear 
combination of atomic orbitals in the Slater-Koster tight- 
binding approximation [sJl- The PFMT method, pre- 
sented in this work, can be classified as a finite difference 
approximation technique. In comparison, an early dis- 
cussion of this method can be found in [Sq , (36j , using a 
basic analytical formulation for the elementary problem 
of the electronic transport across a single chemical defect 
in a model wire nanojunction. 

The combination of PFMT and tight-binding model 
leads to transparent matrix structures which are numer- 
ically solvable in a direct manner, in contrast with the 



iterative matrix algebra characteristic of the NEGF cal- 
culations which require special numerical optimization 
methods [32|- Furthermore, integrated tight-binding ap- 
proximation allows us to discretize the entire model cal- 
culation in the real-space, and equally important, to 
make time economies in the numerical computations. 
Here we note that the tight-binding theory has been used 
successfully in previous work as a basis for electronic 
transport calculations, for various low-dimensional sys- 
tems, such as polymer nanojunctions [i^, carbon nan- 
otubes [4^, and the presently popular graphene systems 
47| . This reinforces our choice of the electronic dynamics 



model 

As an implementation of the PFMT method we con- 
sider in the present work two types of atomic wire nano- 
junction systems at zero bias limit. In the first step our 
model calculations are applied to various lengths of the 
freely suspended monoatomic linear sodium wires (ML- 
NaW), as nanojunctions between one-dimensional semi- 
infinite atomic leads of the same material. These calcu- 
lations are motivated by the fact that the MLNaW nano- 
junctions, which have been investigat ed p reviously using 
other methods 
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12| . [26| . [37| . |38| . are known as 
good benchmark systems for electronic conductance cal- 
culations [2l|, [11]. This initial analysis is hence carried 
out in order to demonstrate the correctness and function- 
ality of the PFMT method by comparing with previous 
results. 

Further, our PFMT method is applied to the diatomic 
linear copper-cobalt wires (DLCuCoW) which are hith- 
erto untreated. The particular interest in these systems 
stems from the fact that they arc mechanically and ther- 
modynamically stable at low temperatures, supported as 
they are on atomically flat Cu(lll) surfaces ^3], [HI, 
[4^ . Our results for this nanojunctions supplement in a 
natural way recent experimental investigations on finite 
Cu-Co wires [1^, by providing a theoretical description 
of the electronic transport properties for these systems. 

The present paper is organized as follows. In section HIl 
we give an insight into the general and essential features 
of the PFMT approach. Section Hill presents the applica- 
tions of our model calculations to MLNaW and DLCu- 
CoW nanojunctions. The numerical analysis yields both 
the transmission and the refiection probabilities for elec- 
trons, as well as the total electronic conductance for the 
considered systems at zero bias limit. The conclusions 
are presented in section ITVl 



II. THEORETICAL MODEL AND METHOD 

The general formulation of the PFMT method is pre- 
sented next. It is convenient to develop it initially, and 
to rewrite it afterwards into the form which is directly 
applied for our model calculations for the MLNaW and 
DLCuCoW nanojunctions. This allows us to introduce 
the PFMT technique in the most detailed and appro- 
priate manner. Furthermore, the general presentation 
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clarifies the mathematical approach which is applied in 
the present work for the system of wire nanojunctions. 
It equally illustrates how the method may be applied in 
general, treating general molecular forms for the nano- 
junction when conjugated between one-, two- or three- 
dimensional leads. The general formulation can also in- 
corporate the interactions between atomic sites beyond 
nearest-neighbors. 

Let us first refer to the Landauer-Biittiker approach for 
the a naly sis of the electronic transport via nanostructures 
[13 J 01 J and divide an arbitrary nano junction system 
into three main parts [s^ , [i^ , namely the left and right 
leads, and the nanojunction domain itself, as is in FigU] 

To model the electronic properties of this system, the 
starting point is to write the general form of the tight- 
binding secular equation, for the total M number of 
atoms per unit cell, and the total L number of atomic 
orbitals per site, as follows 



H.3,.2 H.2^.i H.]_„ H„_i 
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(1) 



Note that we assume spin degeneracy in Eq. ([T]) . In what 
follows, the M by M overlap sub-matrix S„.„/ for a given 

n-th unit cell, incorporates the s[^' ^, overlap integrals 
between the I and V atomic orbitals centered at the sites 
corresponding to the real space vectors and Vm' , re- 
spectively. The unit cell may contain more than one in- 
equivalent site, denoted here by to. Furthermore, the 
corresponding diagonal Hamitonian block matrix H„.„' 
for the n-th unit cell, is composed of binding energy pa- 
rameters eL, and the interaction integrals /ijl' /. Then 

lit ^ liti ^ iiL 

again, the off-diagonal H„^„/ sub-matrix contains only 
1 1' 

m' tsrms which describe interactions between differ- 
ent unit cells. Finally, the c(r„,k) denotes vector of the 
c{yi —Ym — r„ , k) wave function coefScients which corre- 
spond to the / atomic orbital at the lattice site to in the 
ri-th unit cell, k is the appropriate wave vector of the 
Bloch like wave function of the electronic excitation. 

In order to calculate the electronic transport properties 
in the framework of the PFMT method for the nanoscale 
system presented in Fig[l] we have to analyze the elec- 
tronic scattering at the embedded nanostructure. To that 
end, we consider the nanojunction irreducible boundary 
domain, and the interactions of its electronic states with 
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FIG. 1: Schematic planar projection in the direction of prop- 
agation of an arbitrary tliree-dimensional Landauer-Biittiker- 
type system. Tlie unit cells of the nanojunction and two semi- 
infinite leads are denoted by red and blue colors, respectively 
[4^ 1. The corresponding Hamiltonian on- and off-diagonal 
sub-matrices are also depicted. Note that in order to keep 
the figure simple we show the Hamiltonian sub-matrices for 
interactions between different unit cells only for n < n' . 



the electronic excitations incident from the leads to the 
left and right of the domain. The wave function coeffi- 
cients of consecutive unit cells in the leads are defined in 
the sense of the Bloch-Flouqet theorem by the following 
phase relation 



c(r„,k) = zc(r„_i,k). 



(2) 



where z = exp(±ikr„) denotes the generalized Bloch 
phase factor. In our notation z = exp(+ikr„) describes 
the electronic wave function propagating to the right, and 
z = exp(— zkr„) to the left. We should emphasize that in 
the case of multi-channel scattering processes one has to 
consider not only propagating but also evanescent waves 
[33! . The evanescent solutions can be obtained using dif- 
ferent procedures, however an elegant method presented 
previously for phonon and magnon problems [33| shall be 
also applied to the electronic problem. 

For an electron incident along the leads at a given en- 
ergy E and wave vector k, where E ~ E^O^) denotes 
the dispersion curves for the 7 = 1, 2,..., propagating 
waves in the Brillouin zone, the scattering at the bound- 
ary yields coherent reflected and transmitted fields. To 
calculate the Landauer-Biittiker electronic conductance 
via a given nanojunction, we write the linearized equa- 
tions of motion for the consecutive atomic sites within the 
irreducible nanojunction domain, on the basis of Eq.(IT]). 
This procedure generates a set of iV -I- 2M equations of 
motion with N + AM unknown variables, namely the 
wave function coefficients, extending to outside the do- 
main and inside the leads. Since the number of unknown 
variables is always greater than the number of equations, 
such a set of equations cannot be solved directly. 

To resolve this problem, we use the phase field match- 
ing method [ll], [s^, [s^. This implies that, for unit cells 
in the leads distant from the nanojunction boundary, the 
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coefficients may be expressed in terms of an appropriate 
superposition of the eigenstates of the perfect leads at the 
same energy and wave vector. Hence, the wave function 
coefficients for the region outside the irreducible nano- 
junction domain may be expressed as 

c(r„,k) = c^(k)z-" 

r 

+ ^CY{k)z"rj^Y for n ^ -1, (3) 

7' 

r 

c(r„,k) = ^c^'(k)z"f^,-y' for n^N. (4) 

7' 

In equation ([3|) and (jl]), the r-y^ji and i^,^' factors de- 
note the reflection and transmission coefficients for the 
incident electronic wave function corresponding to the 
7 channel, which is reflected or transmitted into the 7' 
channel. Here, by channel we understand one of the total 
number F of the allowed eigenstate solutions for the gen- 
eralized Bloch phase factor in Eq.(l2]), which determine 
the dispersion branches and the evanescent waves for the 
electronic structure of the lead system. The c^(k) term 
denotes hence the lead Hamiltonian eigenvector which 
corresponds to the 7 electronic eigenstate, and serves as 
a weighting factor in equations ([3]) and Q. For a com- 
plete description of the scattering processes it is essential 
to know the evanescent as well as the propagating elec- 
tronic eigenstates fssj . 

Such formulation of the problem is similar to those 
presented in [2^ and [i^. However, in contradiction to 
the mentioned techniques, we use equations ^ and Q in 
order to directly rewrite the equations of motion for the 
region outside the nanojunction domain [s^l, [1^, [36l |. 
In this manner we create an irreducible vector of the wave 
function coefficients for the nanojunction domain and its 
matching regions in the input and output leads, in the 
following way 

6(k) = [z(k) X c"'^"°(k)] + I(k), (5) 

where z(k) is the N + AM by -I- 2M matrix composed 
of the phase factors which are mapped onto the atoms in 
the matching region outside the irreducible nanojunction 
domain. The vector c"°"°(k) incorporates the vector of 
wave function coefficients on the nanojunction domain, 
and the reflection and transmission coefficients for the 
wave function in the matching region in the input and 
output leads. The number of reflection and transmission 
coefficients is equal to the number of scattering chan- 
nels in an appropriately constructed Hilbert space. The 
vector I(k) decouples terms which describe the incoming 
wave. On the basis of Eq.(l5]), we reduce next the resul- 
tant N + 2M set of equations for the irreducible region, 
to derive the following square matrix of inhomogenous 
linearized equations 

[ES - M{z, k)] X c"""'(k) = -i(k). (6) 



S, the overlap matrix, and M, composed of the Hamilto- 
nian elements and appropriate phase factors, are N + 2M 
by -I- 2M square matrices, whereas I(k) is the I(k) 
vector supplemented by the corresponding Hamiltonian 
terms. We observe that, contrary to the Green's function 
methods, Eq.(|6]) allows us to solve the scattering problem 
directly for real energies. 

Moreover, Eq.® yields directly the t^^^i — t^^^i{E) 
and r^.'y' = r^^^i{E) scattering amplitudes, and the wave 
function coefficients in the irreducible boundary domain, 
as a function of the electronic energy E or wave vec- 
tor k. Here E = ^'^(k) denotes the dispersion relation 
for the propagating waves of the electrons in the first 
Brillouin zone of the leads. It is these E = E^{k) elec- 
trons, incident from the perfect leads, which scatter at 
the nanojunction. 

In what follows, the t^.^' — t^^^i{E) and ^'7,7' — 
r~f^^i[E) amplitudes yield the required Landauer-Biittiker 
total transmission and refiection probabilities on the 
nanojunction, T{E) and R{E) 

R{E)^J2^\^'rn'm\ (7) 

7,7' 

TiE)^Y.^Ky(^)\'^ (8) 

7,7' 

summing over all incident, reflected, and transmitted 
waves. To ensure the unitarity character of the scatter- 
ing matrix, the amplitudes are normalized with respect 
to their electronic group velocities, v-y and Vji, for elec- 
trons incident in channel 7 and reflected or transmitted 
in channel 7'. The unitarity condition, T{E) + R{E) = 1, 
corresponds to the necessary conservation of energy. 

Finally, assuming the zero bias limit and spin degen- 
eracy, the overall conductance is then simply written as 

GiN,EF)^GoT{EF), (9) 

where Go = 2e^/h is the conductance quantum. The 
electronic conductance is given at the Fermi energy since 
the electrons at this level give the only significant contri- 
bution to the overall conductance from the Fermi-Dirac 
distribution. 



III. NUMERICAL RESULTS AND DISCUSSION 

In this section, we present the dynamic equations for 
the MLNaW and DLCuCoW nanojunctions on the ba- 
sis of the theoretical discussion presented in section [Hi 
The schematic representations of the freely suspended 
MLNaW and the supported DLCuCoW nanojunctions 
between onc-dimensional atomic electric leads are pre- 
sented in Figl2jA) and Figl2l^B), respectively. For both 
systems each of their atoms is characterized by a one va- 
lence s-type electronic state, and the interactions with 
adjacent atoms have a nearest-neighbor range, which is 



5 



usually proposed for the Na atoms [44|. As for the con- 
sidered Cu-Co systems, it has also been demonstrated, 
in view of previous experimental results [i^l, that their 
electronic dynamics can be described within the s-like 
monovalence approximation, suggesting that one chan- 
nel is sufficient to model the electronic transmission. 

For the numerical calculations of the electronic trans- 
port across MLNaW and DLCuCoW nanojuntions pre- 
sented in this section, the overall Hamiltonian of the 
wire nanojunctions connected to the semi-infinite one- 
dimensional leads can be written on the basis of Eq.([T]) 
in the following reduced form 



H 



n' ,n+l + Sn',n-l), 



(10) 



where we neglect the summation over different electronic 
states per site, and index m is simply replaced by n due 
to the fact that only one atom per unit cell is assumed 
in our calculations. 

In the case of the MLNaW systems (see FigUJA)), 
the Hamitonian elements are assumed to be consistent 
with the Harrison theory convention [i^. In the pres- 
ence of local charge neutrality, all the binding energies 
which characterize the atoms on the MLNaW nanojuc- 
tions are equal to each other, whereas the corresponding 
Ki,n' parameters are distance dependent and defined as 



hn n' 



re 



V- 



(11) 



where 77 = —1.32 is the dimensionless Harrison's co- 
efficient for the s-type orbitals, rUe denotes the elec- 
tron vacuum mass, and d depicts the nearest-neighbor 
interatomic distance. For the MLNaW nanojunctions, 
Eq. (|TT1) allows us to model the effective nearest- neighbor 
interactions by considering different interatomic spac- 
ings. On the other hand, the tight-binding parameters 
for DLCuCoW nanojunctions (see Figl5](B)) are assumed 
after [i^, to be distance independent. 

In particular, for the MLNaW, and DLCuCoW sys- 
tems Eq. (jni can be reduced to the following explicit form 



h-i-2 -^-1.0 
-/io^_iz^ E - '■■ 
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The total transmission and reflection probabilities, for 
the left and right leads, which arc identical atomic wires 
presenting equal group velocities dE/dk for the electrons. 
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FIG. 2: Schematic representations of the: (A) Ai'-atomic lin- 
ear nanojunction made of sodium atoms (light green color) 
between sodium leads (dark green color); (B) A*'-atomic lin- 
ear nanojunction made of pairs of copper (orange color) and 
cobalt (blue color) atoms between copper leads. The cor- 
responding binding energies En for sites n, and the nearest- 
neighbor couplings hn^n' between sites n and n' are depicted. 
Note that in order to keep the figures simple we show the 
nearest-neighbor coupling terms only for n < n' . Addition- 
ally, in (A) three dilTerent interatomic spacings are considered: 
dww for the nanojunction wire, dLL in the leads, and c?wl at 
the contact. The detailed description of the electronic and 
structural properties is presented in section Hill 



are related to the single channel amplitudes t{E) and 
r{E) in the following way 



T{E) = \t{E)f and RiE) = \r{E)\'^ 



(13) 



At this point we would like to draw attention to the fact 
that the PFMT formulation of the scattering problem 
presented in equations ([Slini), allows us to introduce addi- 
tional simplifications to the resultant equations. In par- 
ticular, the boundary terms of matrix [EI — M(2, k)] in 
Eq. ([T2|) . known as a self-energies in the Green's func- 
tion formalism, are real and energy independent in the 
PFMT method. This simplification can be also applied in 
the general case of Eq.([6|) when the interactions between 
neighbor atoms in the leads arc symmetric. 

Note, furthermore, that since the s-state electrons of 
the Na or Cu one-dimensional semi-infinite leads popu- 
late up to half of the available band states, one can show 
by symmetry that Ep is equal to the binding energy of 
the lead atoms. 
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A. Monatomic Sodium wire nanojunctions 

The PFMT model calculation is next applied nu- 
merically for the MLNaW wire nanojunctions between 
sodium leads as in Figj^JA). Our numerical results are 
determined for wires of various lengths, in particular they 
are presented typically here for N € [1,6]. 

As stated in section HH binding energies of Na atoms 
are constant and equal -4.96 eV, while the distance de- 
pendent coupling terms are modeled using Ea.([TT|). The 
interatomic distance on the Na leads Jll = 3.64A is con- 
sidered to correspond to the bulk Na nearest-neighbor 
distance. However, d of Eq. lfTTj) on the MLNaW nano- 
junction takes on a value of dww = 3.39A which is the 
average interatomic distance for the finite Na wires [l^l . 
In contrast, the lead-nanojunction spacing is considered 
for two numerical values, dwL which is equal to 3. 99 A 
and 3.29 A. These numbers correspond to a choice range 
of approximately 20% weaker (WC) and stronger (SC) 
lead-nanojunction couplings, respectively, compared to 
the intrinsic lead coupling. This is a physically reason- 
able range to test. Taking into account the above values 
of the interatomic distances, we note that the N=l case 
is equivalent to the situation of the infinite chain with 
the two bond defects, however this is not true for the 
remaining wire lengths due to the fact that ^ll c^ww- 

In FiglJl the overall transmission T{E) and reflection 
R{E) probabilities are presented as a function of the inci- 
dent electronic energy E, for various wire lengths N , and 
for the two different lead-chain coupling values SC and 
WC defined above. For each sub-figure of FigOl T{E) 
and R{E) are calculated independently, and the unitar- 
ity requirement for the scattering processes is checked by 
the sum T{E) + R{E). Note that Ep is set as a zero 
energy reference in these figures. 

The positions of the transmission spectral resonances 
vary as a function of the lead-nanojunction couplings, as 
expected and as may be seen from Fig|31 Due to the mir- 
ror symmetry and assumed local charge neutrality in the 
considered MLNaW systems, all energy levels, excluding 
the central peak for odd iV, are distributed symmetri- 
cally with respect to the referential Ep, for both even and 
odd N. Analyzing the transmission spectra, the number 
of the resonances is always equal to the number of the 
atoms in the considered MLNaW wire up to = 3. This 
equality is also conserved for A = 4 for the strong SC 
couplings, but reduced to A — 2 for the weak WC cou- 
plings. The A — 2 states are also observed for the five 
odd- and the six cvcn-numbcrcd wires for both SC and 
WC cases. The displacement of the resonances to outside 
the transmission spectra for any A is directly related to 
the lead-nanojunction couplings for the considered sys- 
tem. This displacement happens when the resonance 
goes to a localized state outside the leads monovalent 
s-type conduction band. This effect is equally observed 
for resonances and localized vibration states on atomic 
nanojunctions [11]. Further, when the AILNaW wire in- 
creases in length this displacement is quicker for the WC 
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even-number. Ep is set as a zero energy reference, and the 
unitarity condition is represented by the sum S{E). 
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connected wires than for the SC ones. 

As mentioned previously, the only significant contri- 
bution to the overall conductance comes from electronic 
states in the neighborhood of the Fermi level Ep. It 
is hence noticeable that for the odd-number wires, one 
of the resonant maxima appears systematically at Ep, 
which yields as a consequence a conductance maximum 
for the corresponding systems. In contrast, the even- 
numbered wires do not present any maxima at Ep, which 
outcome results in successive minima for the conductance 
of these systems. 

The total conductance G{N,Ep) for the considered 
MLNaW wires is presented in FiglU as a function of 
the wire length iV, reflecting the maxima and minima 
of the transmission spectra at Ep, and exhibiting the 
usual conductance oscillations, here with a two-atom pe- 
riod. This situation relates directly to the interaction of 
the free electrons in the lead band with the electronic 
states of the nanojunction wire. In particular, since the 
a state for a sodium atom can take two electrons by the 
Pauli principle, one notes that there are half-empty states 
available for the odd-number wires while non available 
for the even-number wires. Since the Fermi energy Ep 
is equal to the binding energy of the Na atoms, the odd- 
number wires present a greater local density of states 
LDOS accessible at the Fermi level D{Ep), compared 
to that for the even-number wires. This ensures that the 
hopping mobility is unimpeded for the odd-number wires, 
whereas this is not the case for the even- number wires. 
The corresponding electronic transport is consequently 
a maximum at ^ Go for the odd-number MLNaW wire 
nanojunctions, for both the weak WC and the strong 
SC couplings. The electronic conductance decreases for 
the even-number wires, and this decrease varies with the 
lead-nanojunction couplings, it is ~ O.OGSGo for SC and 
^ O.22IG0 for WC, respectively. Furthermore, one ob- 
serves that the conductance of the even- number MLNaW 
wire nanojunctions increases slightly but monotonically 
with the increase of the wire length. This situation may 
be explained qualitatively by arguing that the electric 
leads have a bigger influence on the LDOS at the Fermi 
level for the short wires than for the longer ones. An in- 
teresting discussion on the problem of how the different 
types of leads may influence the overall electronic con- 
ductance can be found in [l^ . 

In summary, the electronic conductance determined in 
this subsection presents hitherto unknown detailed trans- 
mission spectra, and confirms the well known even-odd 
conductance oscillation effect [soj , which is previously re- 
ported for the MLNaW nanojunctions @, as well as for 
other wire systems [l^l, [16j . We also note that our 
results are in agreement with those obtained on the basis 
of the first-principle calculations in other model calcula- 
, [33] ■ This clearly confirms the validity of 
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the electric leads. 



B. Diatomic Copper-Cobalt wire nanojunctions 

The second considered system in this work concerns 
the supported DLCuCoW wire nanojunctions presented 
in section II. Finite - Cu - Co - wire systems have been 
experimentally prepared previously as grouped adatoms 
on the fiat Cu(lll) surface by Lagoute et al. [i^l, with 
the use of the low-temperature scanning tunneling mi- 
croscope technique. Initially, these systems have been 
developed as magnetic/nonmagnetic finite atomic wires 
for the investigation of their potential magnetic proper- 



Hi, 



tions 

the PFMT approach. The minor quantitative differences 
between our results for the conductance and those pre- 
viously reported, for the same A'^-number MLNaW wire 
systems, appear mainly due to the assumed structures of 



ties. Subjects concerning their spin dynamics |5l| . their 
storage potential [s^l, and the observation of the Kondo 
effect in their structure (ssj , have been indeed addressed 
recently. However, their electronic transport properties 
have not been considered previously, despite their intrin- 
sic interest under a DLCuCoW configuration, and their 
mechanical and thermodynamical stability at low tem- 
peratures. Note in this respect the increasing interest 
in the stability of such supported wire systems [U for 
arbitrary ranges of temperature. 

As it was stated in section II, the electron dynam- 
ics of the supported DLCuCoW nanojunctions can be 
discussed within the s-like single band effective tight- 
binding model proposed in (40| . This possibility arises 
from the existence of unoccupied confined quantum 
states formed along the Cu-Co wire nanojunctions by 
linking the orbitals of spz character. Such confined states 
have been observed in other supported wire systems like 
the monatomic Cu [4l[ and Au js!] , the diatomic Au-Pd 
nanojunction systems [H^, and Fe nano-islands j57| . 

In this context the possibility of using semiconductor 
surfaces which present a band gap is equally interesting 
as a support for the atomic wires. These surfaces may 
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FIG. 5: (A) - (D): The selected spectra of the transmission T{E) and reflection R{E) probabilities for respectively 2, 7, 12 
and 17 copper-cobalt atomic pairs in the DLCuCoW scattering region. The unitarity condition is represented by the sum 
S{E). (E): Surface plot of the transmission probabilities as a function of energy and number of the Cu-Co atomic pairs in the 
DLCuCoW scattering region. For all subfigures the Ef is set at zero of energy. 



electronically decouple from the supported wire orbitals 
in certain cases, allowing one to study the electronics of 
stable one-dimensional wire systems (5J]. 

From our point of view the diatomic nanojunction 
wires may be used to control the value of the overall 
electronic conductance in low-dimensional nanoelectronic 
circuits. In particular, the foreign cobalt atoms injected 
directly into the monatomic copper wire, act as control- 
lable chemical defects which modify the initial system 
and its electronic properties as a consequence. Further- 
more, the substrate supported configurations of the lead 
- DLCuCoW wire nanojunctions - lead, suggest a greater 
mechanical and thermodynamic stability than for the 
freely suspended wires, which is one of the most impor- 
tant features required for future nanoelectronic devices. 

Our model calculations for the electronic quantum con- 
ductance of the DLCuCoW systems are presented for 
wires composed of different numbers of periodically ar- 
ranged Cu-Co atomic pairs, as in FiglJl^B). The length 
of the DLCuCoW wire may characterized by the number 
-^CuCo of such pairs, and the results in this subsection 
are hence presented in terms of this latter variable for a 
given DLCuCoW nanojunction. 

In the case of the DLCuCoW nanojunctions, we as- 
sume the Hamiltonian elements of Ea. ([TU)) to have the 
values proposed in [i^l . Following these studies, the bind- 
ing energy 3.31 eV is considered to be the same for all 
the copper atoms, whereas for cobalt it is 2.96 eV. The 
coupling terms between copper nearest-neighbors on the 
leads, and between copper and cobalt nearest-neighbors 
on the wire nanojunction, are close to each other with 
values of —0.95 eV and —0.94 eV, respectively. 

It is important to note that the iVcuCo = 1 nanojunc- 



tion wire is a very particular case since it has quite a 
different symmetry from all other lengths and is open to 
different structural definitions; we have calculated the to- 
tal conductance for this particular case to be G.OGSSGo, 
but will not use it in our discussion of the properties of 
the conductance G(A^cuCo, Ep) for a general wire length. 

In Figl5jA)-(D) we present the detailed transmission 
T{E) and refiection R{E) spectra of the DLCuCoW 
nanojunction, for the selected values of iVcuCo = 2, 7, 
12, and 17. The resonance maxima and minima can be 
observed as for the MLNaW wires, where Ep is again set 
as a zero energy reference. However, the transmission 
and reflection spectra are not symmetric with respect 
to this reference, due the slight difference between the 
cobalt and copper binding energies. The most impor- 
tant conclusion concerns the value of the transmission 
probabilities close to Ep. The transmission spectra in 
FiglSIA)-(D) show strong resonant minima for both the 
even- and the odd-number wires. In particular, we ob- 
serve that the value of the transmission at the Fermi level 
decreases when the length of the effective nanojunction 
increases. This trend is summarized in FiglSl^E) for the 
ensemble of the transmission spectra for A^cuCo values 
ranging from 2 to 17. At the upper limit the conduc- 
tance goes to G{NcuCo = ^7,Ep) ^ IO^^Gq which is 
sufficient for the purpose of the present calculation. 

The overall electronic conductance G{NcuCo, Ep) as a 
function of the DLCuCoW wire length in terms of Cu-Co 
atomic pairs, is presented in FiglHl along with a fitting 
exponential function 

GiNcnCo.Ep) = Gi2,Ep)cxp[~aiNcnCo - 2)] 
for iVcuCo > 2. (14) 
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Fig |6] resumes the transmission spectra at the Fermi en- 
ergy, it shows the strong monotonic decay of the elec- 
tronic conductance with the increase of the number of 
Cu-Co pairs of the nanojunction. In the corresponding 
Eq.dUl), G{2,Ef) = 0.8760Go is the conductance of a 
wire made up of two Cu-Co pairs, and a ~ 0.1630 denotes 
the decay constant. Eq. p^ allows one hence to read the 
numerical values of the conductance G(A^cuCo, Ep) of the 
DLCuCoW wires for iVcuCo > 2 using a simple formula 
which is useful in the analysis of potential applications of 
corresponding nanocircuits. This behavior is fundamen- 
tally different from the typical conductance oscillations 
observed for odd- and even-number monoatomic wires. 



> 
>< 

CD 

c 
LU 









- 








A=0.35 eV - - 


















0.0 I 

2 4 6 8 10 12 14 16 18 

Number of Cu-Co pairs 



FIG. 6: The total electronic conductance G{NcuCo, Ef) in 
units of Go, read at the Fermi level Ef as a function of the 
number of Cu-Co pairs on the DLCuCoW nanojunction wire. 
The PFMT results are represented by open circles, and the 
fitting function by the red curve. 



To understand the decrease of the conductance in 
FigEl one has to take a closer look at the band struc- 
ture of the infinite -Cu-Co- diatomic chain which is the 
infinite limit of the DLCuCoW nanojunction. We have 
calculated the band structure for such an infinite di- 
atomic wire as in Fig 171 This presents a band gap of 
Z\=0.35 eV and renders the system theoretically insulat- 
ing. This band gap, which corresponds directly to the 
difference of the binding energies for cobalt and copper 
atoms, presents a differential ~ 10^^ with respect to the 
Cu binding energy. It is an effective potential barrier 
for electrons between successive Cu and Co sites, and 
results in a small but cumulative decrease of the conduc- 
tance as the DLCuCoW wire nanojunction increases in 
length, when comparing to the theoretical maximum of 
Go for the pure atomic Cu wire. This direct interpre- 
tation is confirmed by the exponential form of Eg. ([Ml), 
where exp[— a(A^ — 2)] = exp(— a)^~^, expressing the 
monotonic decay of the electronic conductance with the 
increasing length of the DLCuCoW nanojunction. 



FIG. 7: Band structure of the infinite variant of the DL- 
CuCoW wire nanojunction over the first Brillouin zone. The 
band gap (Z\=0.35 eV) is marked by two dotted lines, between 
the copper (ecu) and the cobalt (eco) binding energies. The 
Fermi energy for the half-filled bands is represented by the 
dashed line. 



IV. SUMMARY 

In this work we present a general method based on 
the algebraic phase field matching theory (PFMT), to 
calculate the electronic quantum transport across nano- 
junctions of arbitrary form between material leads of n- 
dimensions, n = 1, 2 or 3. This PFMT method is ap- 
plied in particular to calculate the quantum electronic 
conductance across wire nanojunctions held between one- 
dimensional electric leads. The presented PFMT formal- 
ism provides a compact and efficient approach for the 
analysis of the electronic quantum transport for a wide 
range of wire nanojunctions (ssj . 

The model calculations are carried out for nanojunc- 
tions made up of sodium atomic wires (MLNaW) between 
leads of the same element. Our calculated results for the 
electronic quantum transport are shown to be a function 
of the physical parameters of the wires, and exhibit the 
well known oscillation behavior for even- and odd-number 
wire lenghts. They are also in agreement with the first- 
principle results [lH , , [111 . This clearly confirms the 
validity of the presented PFMT model. 

Our model calculations are also carried out for the di- 
atomic copper-cobalt wires (DLCuCoW) held between 
copper leads. This system is selected for its mechani- 
cal and thcrmodynamical stability at low temperatures, 
which is one of the important features required for future 
nanoelectronic devices. In contrast to the MLNaW sys- 
tem, the electronic quantum transport of the DLCuCoW 
wire nanojunction, at the Fermi energy, exhibits expo- 
nential decay with increasing wire length. This behavior 
is explained notably in terms of the band gap of the infi- 
nite variant of the DLCuCoW nanojunction wire. A re- 
lation, depicting analytically the numerical values of the 
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conductance G{NcuCo, Ep) of the DLCuCoW wires for 
any length A^cuCo > 2, is given and may help to simplify 
the analysis of potential applications for corresponding 
nanocircuits. 

The presented PFMT approach, a transparent and 
time-saving method, should be especially interesting for 
the treatment of complex systems presenting multichan- 
nel conductance, and for the treatment of the elec- 
tronic quantum transport of nanojunctions which can ex- 
hibit electron-electron and electron-phonon interactions 
at high temperatures or due to structural disorder. 
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